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Abstract 

We consider non- interacting particles subject to a fixed external potential V and a 
self-generated magnetic field B. The total energy includes the field energy /3 J and 
we minimize over all particle states and magnetic fields. In the case of spin-1/2 particles 
this minimization leads to the coupled Maxwell-Pauli system. The parameter (3 tunes 
the coupling strength between the field and the particles and it effectively determines 
the strength of the field. We investigate the stability and the semiclassical asymptotics, 
/i — 7- 0, of the total ground state energy E{/3,h,V). The relevant parameter measuring 
the field strength in the semiclassical limit is k = (3h. We are not able to give the 
exact leading order semiclassical asymptotics uniformly in k or even for fixed k. We do 
however give upper and lower bounds on E with almost matching dependence on k. In 
the simultaneous limit /i — ?• and k — ?• oo we show that the standard non-magnetic Weyl 
asymptotics holds. The same result also holds for the spinless case, i.e. where the Pauli 
operator is replaced by the Schrodinger operator. 
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1 Introduction 

An important problem in spectral analysis is to determine or bound the sum of the negative 
eigenvalues of a Schrodinger operator —A — V{x), i.e.. 



under appropriate conditions on the potential V . We use the convention that x_ = (x)_ = 
min{x, 0} when x is either a real number or a self-adjoint operator. This problem is of 
particular interest in quantum mechanics as it gives the ground state energy of a gas of free 
fermions moving in the exterior potential V . 

A generalization of this problem is to consider not only a potential V but also an exterior 
magnetic field given by the vector potential A. The corresponding magnetic Schrodinger 
operator is (— zV + A)"^ — V{x) acting in L^(M'^). A further generalization is to consider the 
particles as having spin— ^ and introduce the magnetic Pauli operator [cr ■ (— iV + Ayf' — V, 
where in d = 3 dimensions cr = (cti, (T2, da) denotes the vector of 2 x 2 Pauli matrices. The 
Pauli operator acts in L^(M^;C^). Much work has gone into understanding the semiclassical 
asymptotics of the sum of negative eigenvalues, i.e., the asymptotics for small h > of 



It is well known that under appropriate conditions on A and V the leading behavior as h 
tends to zero is given by the Weyl formulas 



respectively. Here the factor of 2 on the second integral is due to the spin degrees of freedom 
in the Pauli operator, i.e., the fact that it is a 2 x 2- matrix valued operator. Note that the 
limiting semiclassical behavior is non-magnetic, i.e. fixed magnetic fields do not influence 
the leading order semiclassics. For simplicity we will consider the d = 3 dimensional case 
only and we denote the Schrodinger operator T^{A) = {—ihV + A^ and the Pauli operator 
Tf{A) = [cr ■ {-ihV + A)]^. The magnetic field is 5 = V x A. 

In this paper we will address a related and equally important issue, namely the case when 
the magnetic field is not a fixed external field, but the self-generated classical magnetic field 
generated by the particles themselves. 



Tr (-A - V{x)) 



Tr {{-ihW + A)^ -V{x)) 



or Tr ([cr- (-z/iV + - 
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We will consider the external potential y to be a fixed (given) function in M^, we assume 
V e Ll^^{M.^) and we will always work with the flat Euclidean metric. The vector potential 
A will be optimized to minimize the total energy consisting of the energy of the particles and 
the field energy 

j B^ = j |Vx A|2 

(we use the convention that unspecified integrals are always on w.r.t. the Lebesgue mea- 
sure). The problem we consider is thus to determine the energy 

E^'^(p,h, V) =inf 

A 

for /3,h > 0, where the infimum runs over all vector fields A e H\R^; R^); in fact mmimizmg 
only for all A G C^(M'^;M'^) gives the same infimum. See Appendix Rl for a discussion of 
equivalent variational spaces for this energy and for the precise definition of the operator 
T^'^(A) — V and the sum of its negative eigenvalues. We will omit the superscripts S,P, when 
making general statements valid for both the Schrodinger and Pauli cases. 

Here (3 is an additional parameter setting the strength of the coupling of the particles to 
the field. In a given physical system the values of h and /3 are given, but as is standard in 
semiclassical analysis we leave them as free parameters. Formally /3 = oo corresponds to the 
non-magnetic case; smaller /3 means that a larger effect of the magnetic field is expected. 

The Euler-Lagrange equation corresponding to the variational problem ([T]) above is 

PVxB = Ja, (2) 

where Ja is the current of the Fermi gas, which in the Schrodinger case is 

Ja{x) = -Re [(-zW + A)l(_oo,o](T;,(A) - V)] {x,x) 

and in the Pauli case is 

Ja{x) = -Re [Tr {(t{ct ■ (-zW + A))l^_^,o]iniA) - V))] (x, x). 

In other words the Euler-Lagrange equations are the non-linear coupled (time independent) 
Maxwell-Schrodinger or Maxwell-Pauli equations where we deliberately ignored Gauss' equa- 
tion for V in order to obtain a general result. In the application to large atoms discussed 
below we will consider the special case when V solves Gauss' equation. 

There are three natural questions about the energy E{l3,h,V). First of all we may ask 
whether the energy is finite, i.e., not negative infinity. We refer to this as stability. In Theo- 
rem [2TT] we give bounds on the energy in the Pauli case under essentially sharp assumptions on 
V. The corresponding results for the Schrodinger case are well known and are also discussed 
in Section [21 



TT{T^f{A)-Vy 



(3 



IV X A\ 
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The second natural question is whether the inclusion of the self-generated magnetic field 
will actually lower the energy at all. For both the Pauli and the Schrodinger cases this is indeed 
true (see Appendix [B]), i.e., there exist potentials V and parameters /3, h (not necessarily the 
same for the Pauli and Schrodinger cases) such that 

E{(3,h,V) < E{(3 = oo,h,V). (3) 

The third and, actually, main question we address is how the inclusion of the magnetic 
field influences the semiclassical asymptotics {h — )■ 0) for the sum of the negative eigenvalues. 
Standard semiclassical results typically assume that the physical data (external potential, 
vector potential, metric etc.) are smooth. If as above these data arise as self-generated and 
thus determined internally via a variational principle the smoothness is not a-priori given. 

In Theorem 12. 21 we establish under the appropriate assumptions on V that the semiclassical, 
i.e., h ^ 0, asymptotics of -E(/3, h, V) in the case when — )■ oo simultaneously with /i — )■ 0, is 
given by the standard non-magnetic Weyl formula. The case of large Ph is the case of greatest 
physical interest as the magnetic field in general gives rise to a small effect. In |EFS2] we give 
improved estimates on the error term to the non-magnetic Weyl term. 

It is however at least of mathematical interest to understand the behavior of E{l3,h,V) 
also when ph is smaller, i.e., when the effect of the magnetic field is greater. Unfortunately, we 
have not been able to establish the exact semiclassical asymptotics in this case. We do show 
however in Theorem 12.11 again under appropriate assumptions on V that the Pauli energy 
h^E^[P,h,V) is bounded from above and below by functions of k := Ph. These functions 
have almost matching asymptotics in the regime of small ph. In proving the upper bound we 
rely heavily on the construction of zero- modes in |ESlj . The similar result for the Schrodinger 
case is well-known and, in fact, here the bounds do not depend on n (see ([6]) below). 

Theorem 12.11 moreover shows that for sufficiently small k, = ph (depending on V) the 
magnetic and non-magnetic energies are different, i.e., holds, since the former is inde- 
pendent of K, while the latter scales at least as k~^~^'^ for any e > 0. In summary, we have 
established that P ~ is the correct threshold for the Pauli operator to observe the effect 
of the magnetic field in the leading order semiclassics. 

Our results will be stated and proved for the Pauli operator and we will remark the 
modifications for the Schrodinger case. 

1.1 Applications to large atoms 

One of the main applications of precise semiclassical estimates is to investigate the ground 
state energy of large atoms and molecules. It is a celebrated result of Lieb and Simon |LSj 
(see also [L]) that the energy of a neutral atom or molecule with nuclear charge Z obeys the 
Thomas-Fermi asymptotics, — {const. )Z'^^^ , in the large Z limit. The subleading correction to 
order Z^ is known as the Scott correction and was established for atoms in [HI ISWlj and for 
molecules in [IS] (see also |SW2l ISW3[ ISSj). The next term in the expansion of order Z^^^ 
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was rigorously established for atoms in |FSj . The large Z asymptotics can be viewed as a 
semiclassical limit with h = Z~^^^ being the semiclassical parameter. The Scott term thus 
corresponds to the next order semiclassical estimate. In the results mentioned so far on the 
large Z asymptotics, magnetic interactions are ignored. 

Magnetic fields in this context were first taken into account as external fields, either a 
homogeneous one ^ ILSYll ILSY2j (see also [Soblj ISob2t ISob3j for improved semiclassical 
estimates and [Ivll IIv2j for inclusion of the Scott correction) or an inhomogeneous one [ESlj 
but subject to certain regularity conditions. Self-generated magnetic fields, obtained from 
Maxwell's equation ([2]) are not known to satisfy these conditions. In |ES3] the validity of 
Thomas-Fermi theory was extended by allowing a self-generated magnetic field that interacts 
with the electrons. This means that the focus was on the absolute ground state of the system, 
after minimizing for both the electron wave function and for the magnetic field. Without going 
into details we mention that the Thomas-Fermi theory can be viewed as the semiclassical 
approximation to our Maxwell-Pauli system but with the Gauss' equation included. It was 
shown that the additional magnetic field does not change the leading order Thomas-Fermi 
energy. This holds if Za"^ is sufficiently small, where a is the fine structure constant; for large 
values of Za"^ the system is unstable (see [ES3j for more details). 

The semiclassical problem corresponding to a self-generated magnetic field is exactly of 
the type ([1]) we discuss in the present paper for (3h — )■ oo. In order to establish the Scott 
correction for self-generated magnetic fields, it is necessary to establish the semiclassical ex- 
pansion of £'(/?, h, V) up to subleading order if is bounded from below. Such an improved 
semiclassical estimate is proven in a separate paper |EFS2] and the Scott term asymptotics is 
proved in |EFS3j . 

Acknowledgment: JPS thanks R. Seiringer and C. Hainzl for fruitful discussions. 



2 Results 

In this section we will state and discuss our two main theorems. The proofs are given in the 
following sections. 

Theorem 2.1 (Stability bounds for the Pauli energy). Assuming V G L;^„^(R^) and [\^]+ G 
L^/^(]R^) n ./^"^(M^), we have, for the Pauli operator 

h^E^'iP, h, V)>'Cl [V]f - cm-' I [V]i. (4) 

for C > and all h G (0, oo) and (3 G (0, oo]. On the other hand if V E Cq(M'^) then for all 
< £ < 1/3, /i G (0, oo) and P G (0, oo] 

h^E^iP, h, V) < -C I [V]T - C,m-'^^' [ [V]\-' + Mh) (5) 
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for positive constants C and C^, the latter depending on e, and with an error function h h-> 
£v{h), depending only on e, h,and V and satisfying limh^Q £y{h) = 0. 

For reference, we mention the analogous result for the Schrodinger case, 
we have for all /i > and /3 G [0, oo] 

-cj [V]f < h'E'{/3, h, V) < I [V]f + £y{h) (6) 

with some positive constant C and with an error function depending on V and h and again 
satisfying \im.h^Q£v{h) = 0. The lower bound in (E]) is the classical Lieb-Thirring inequality 
[LTj . which holds also for magnetic Schrodinger operators (see e.g., [S]) . The upper bound 
is achieved by setting A = and using the standard Weyl semiclassical estimate for the 
non-magnetic operator. 

Not only does Theorem 12.11 give the almost sharp asymptotics of h^E^{l3, h, V) in the limit 
of small K = Ph, it also essentially gives the optimal condition on which L^-norms of V are 
needed to bound the energy h^E^{Kh~^, h, V) uniformly in h for small h and fixed k. From 
the lower bound in the theorem the L^/^ and norms bound the energy. Conversely, from 
the upper bound all norms with p G [5/2, 4) are needed since if Vn G Cq is a sequence such 
that — > oo as n -> oo for some p G [5/2,4), there is a sequence hn tending to zero 
such that h\E^{Kh~^, hn, V) — oo. 

Remark 1. We can also ask a slightly different question: What L^'-norm condition on V 
will ensure finiteness of the energy, but not necessarily a hound in terms of this norm? For the 
Schrodinger operator the answer is that V+ G L3/2(r3) 

ensures finiteness of E^{l3,h,V) and 
this is essentially sharp as far as the local singularity is concerned, since for the critical case 
V{x) = c|x|~^ the energy is finite if c < h'^/4 and infinite otherwise. For the Pauli case we do 
not know a similar sharp result. 

If we ask instead for finiteness of only one eigenvalue, i.e., the one-electron energy 



E^'''{/3,h,V) = mf 



A 



inf Spec(T^^'^(A) - V) + (3 J \V x A\^ 



(7) 



we have a sharp result for both Pauh and Schrodinger. Of course EQ{f3,h,V) > E{f3,h,V) 
and hence instability for Eq implies instability for E. The situation for Schrodinger is exactly 
the same for Eq and E^, i.e. if V+ G L^/^(M'^) then Eq is finite and for the critical case 
V{x) = c\x\~^ the energy is finite if c < h'^/4 and infinite otherwise. 

For Pauli Eq is finite if G L^(]R^) and for the critical case V{x) = c\x\~^ there exists 
a critical value 7cr > such that E^ is finite if c < 7cr/3/i^ and infinite if c > 'jcrPh'^- This 
is essentially contained in |FLLj and we shall review it briefly in Appendix O We will also 
discuss in the appendix that the sum of all eigenvalues E^ likewise remains bounded for the 
cutoff Coulomb potential V{x) = [c\x\~^ — 1]+ if c > is small enough. 

Finally, we give the result of the exact spectral asymptotics in the case of weak magnetic 
fields. 
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Theorem 2.2 (Semiclassics for weak fields). Assume that V G L^/^(]R^) fl L'^(R^) then 



lim h^E^{f3, h, V) = lim h^E^{oo, h, V) = ^ [ \V{x)ffdx. (8) 

''^0 h^o IStt / 

Likewise for the Schrddinger case we have for V E L^/^{R^) 

lim h^E^{(3, h, V) = lim h^E^{oo, h, V) = ^— [ [V{x)ffdx. (9) 

I, V ^ ^ 



/3h—^oo 



Remark 2. This result is a strengthening of Theorem 1.3 from |ES3j . where the same 
conclusion was proved under the condition > c > 0. 

Remark 3. We conjecture that in case of the Schrodinger operator the Weyl term is the 
correct asymptotics uniformly in /3, i.e. iiVE L^/^(M'^), we conjecture that 

lim h'E''{/3, h,V) = --^ [ [V{x)]fdx (10) 

holds uniformly in /3 G [0, oo). The upper bound ([5]) shows that the same statement cannot 
hold for the Pauli case. 

Remark 4- Subleading error estimates are established first in |EFS2j and later in |Iv3] . 

Remark 5. The subleading error estimates in |EFS2] are used in [EFS3j to give the two 
term energy asymptotics, i.e., up to the Scott correction, in the large nuclear charge limit for 
atoms and molecules in a self-generated magnetic field. 



3 Stability bounds: Proof of Theorem 12.1 



The lower bound in Theorem 12.11 establishes the stability of the system, i.e., boundedness from 
below of the energy. As the sum of the negative eigenvalues for the Pauli operator T^{A) — V 
goes to minus infinity (if 7^ 0) as the magnetic field increases (e.g. for a constant magnetic 
field |AHSj ). the addition of the field energy is necessary for stability. Moreover, stability will 
also require sufficient decay of the potential V at infinity and control of the local singularities. 

The basic tool for the proof in the Pauli case is the magnetic Lieb-Thirring inequality from 
[LL5] : 

Theorem 3.1 (Magnetic Lieb-Thirring inequality |LLS] ). There exists a universal constant 
C > such that for the Pauli operator Th{A) — V with a potential V G L;^q^(M'^) and V+ G 
L^/2(M3) n L^iR^) and magnetic field B = V x Ae ^^(M^) we have for all h > 

Tr [T,{A)-Vl>-Ch~'J [Vff-C[h~^j \B\^f\ J [Vjlf". (11) 
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Proof of the lower bound in Theorem \2.1\ Using the magnetic Lieb-Thirring inequality 
(ITT]) , we have 

Tr [r,(A)- + j |Vx A|2 

> -Ch-' j [v]f-Ch~'''^(^ j [V]Xf"(^ j B^Y' + P j B' (12) 

> -ch~' j [vf/^ - ch~'{f3hy' j [v]X, 

where we set i? = V x A and optimized over j B^. □ 
3.1 Upper bound in Theorem 12.11 

In order to construct a trial state that will give the upper bound in Theorem 12.11 we first 
use the method of [ES2] to show that there exist compactly supported magnetic fields and 
corresponding Pauli operators with zero-modes with arbitrarily fast decay. The original con- 
struction of zero-modes in |LYj will neither lead to compactly supported magnetic fields nor 
arbitrarily fast decaying modes. 

Proposition 3.2. Given m G N. There exist a smooth magnetic field of compact support 
B : M.' ^ M.' with a corresponding smooth vector potential A : M.' ^ M.', a smooth and 
non-vanishing ip G C°°(]R^; C^), and a constant C > such that 

(T- (-iV + A)^ = 

and 

\i){x)\ < Clxl-^'-K 

Proof. We use the construction and notations from |ES2j . Consider the map 
^■.R' 3x = (xi, X2, xs) ^ 2 ^^^^*-'^ ^ G C U {oo}. 

Xi + 1X2 

This map is $ = r2 o o r3'^|K3, where ra : — )■ U {oo}, and r2 : §^ — )■ U {oo} are 
stereographic projections (we have identified C and M^) and : S'^ — )■ §^ is the Hopf map (see 
Lemma 34 in [ES^ l Note that 

l + i|$(x)p = ii^^^>NVl6. (13) 
4 xi + X2 

Consider the real 2-form on C, which by stereographic projection pulls back to (1/4 times) 
the volume form on i.e., 

u = -{1 + \z\'^/4:)~'^idz A dz. 
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We have J^uj = it. For g a real smooth compactly supported function g G C^(C) we define 
the 2-form 

(3s = ^*{guj) 

on M^. Note that is a closed 2-form. We can therefore define a (divergence free) vector field 
B such that S ■ (X x Y) = /33(X, Y) for all vector fields X and Y. By ([12]), and hence B 
has compact support in M^. Moreover we will assume that (27r)~^ f^gu = m + ^. 

Consider also the real 2-form (32 = {g — l)uj on C. It satisfies (27r)~^ (32 = m. If we set 

h{z) = TT-i [ \n\z- z'\^(32{z'), 
Jc 

then (32 = da2, where 02 is the real 1-form 02 = 2Re (^^dzh(z)dz^ . According to the Aharonov- 
Casher Theorem (See Theorem 37 and Appendix A in |ES2] ) the magnetic Dirac-operator on 
C corresponding to the metric with volume form u and conformal to the standard metric and 
one form (magnetic vector potential) «2 has a (unnormalized) zero- mode, i.e., element in the 

kernel, of the form f{z) , where 

/(2) = (l + |^|V4)V2exp(-M^)/4). 

In the standard metric on C the same would be true with a zero-mode without the prefactor 
(1 -|- |zp/4)^/^, which comes from the conformal factor according to Theorem 23 in |ES2] . 

Note that \f{z)\ < C|2;|^~™. In fact, there will be m zero- modes, but we need only the one 
given above which has the fastest decay. 

Let us turn to the construction of the zero- mode in M.^. Let f2 : — )■ be the vector 
field corresponding to the 2-form $*(a;), i.e., such that O ■ (X x Y) = a;($*(X), $*(Y)). We 
may choose a smooth unit vector field ^ : M'^ — !■ such that cr ■ = |r2|^. In fact, an 
explicit choice is 

e(x) = (i + |xiV4)-^/^ + (J)- 

According to Section 8 in |ES2] we can find a smooth A with V x A = B and such that 
cr ■ (— zV + A)t(! = 0, where 

ij{x) = {i + \x\yfmxmix). 

Here again (1 -|- |xp)~^ is the conformal factor coming from the stereographic projection T3 
(see Theorem 23 in jES2] ). Thus we see that \ip{x)\ < Clxj-'^-^ □ 

We can use this proposition to construct a low energy one-electron state localized in a ball. 

Proposition 3.3. For all < 6 < 1 there is a constant Cs > such that for all h, (3 > and 
all balls Br with radius R we can find a smooth magnetic field B supported in Br such that for 
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all smooth magnetic vector potentials A with V x A = B in we can find an L"^ -normalized 
ijeC^iBnX^) satisfying 



J \(T ■ {-ihV + A)i)\^ + (3 J B^ < Csh^P^-^R'^'^ 



Proof. This is a simple localization and scaling argument based on the result in the previous 
proposition. First note that we may assume that is small enough. In fact, if > c for 
some constant c then h"^ f3^~^ R~^^^ > c^~^h^R~'^ and an upper bound of the form Ch^R ^ can 
be achieved by choosing A = B = Q. 

Without loss of generality we may assume that the ball Br is centered at the origin. Choose 
an integer m such that (2m) ~^ < 5. Let be a zeromode as constructed in Proposition 13. 2[ 
i-e, 

(T- {-iW + A)^p = 

with decay |'?/'(a;)| < Clxj""*"^ and corresponding magnetic field B = V x A of compact 
support. We may assume that tp is normalized. For £ > we define 

(x) = , Ae (x) = r^hA{x/i) . 

Then Bi{x) = hi-^B{x/i) and 

cr ■ (-i/iV + Ae)'ilJi = 0. 

We can assume B to be supported in a ball of radius 1 centered at the origin (otherwise we 
rescale as just explained). Hence Bi is supported in a ball of radius £. 

Choose X ^ C^(]R'^) with support in the unit ball centered at the origin and such that 
< x{x) ^ 1 for 3.11 X and x{^) = 1 if < 1/2. Set 

ip{x) = J\f'~^x{x/R)^e{x). 

1 /2 

where the normalization constant is A/" = (/ \x{x / R)4'£ix)\'^dx) . Then ^ is supported in 
Br. Moreover, 1 - C{l/Rf'^-^ < A^^ < 1, and thus if > 2£ 

j \<T ■ i-zhV + Ae)tp\' = h'AT-' j {y{x{x/R))fm^)\''dx 

< Ch^R-\i/Rf"'-\ 

Hence 



\(T ■ {-ihV + Ae)ij\^ + 13 I Bl< Ch^R-^i/R)^""-' + Ch^pr^ 

= Ch^PR-\PR)-^ 
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with the optimal choice i = CR^j^R)"^ < i?/2 if /3_R is small enough. 

Finally, if A is any smooth vector potential with V x A = 5^ in Br then we can gauge 
transform, i.e., find a smooth </> : Bji — j- M such that Ai = A + V(j) in Bfi. Then 

cr ■ {-ihV + A)e''''''l'^lj = e'^^'^a- ■ {-ihV + v4^)V^ 

and thus the above bound holds with ip replaced by e"*'* ^'^ip. □ 

Proof of the upper bound in Theorem \2.1\ By choosing A = we can always achieve a Weyl 
upper bound 

h^'E^P, h, V)<-C J [Vff + Svih). 
Let us now show that if Wh)-^+^' J[V]X-' > J[Vff then we can achieve the bound 

h'E{f3, h, V) < -C,{f3hy'^'' [ [V]i-' (14) 



for h small enough depending on V and e. Divide space into cubes of side length For 
each cube consider the minimal value Knin of V. Then in the cube V < Vmm + v^37i|| W||oo- 
If Vmin < V^||VV||oo "we do nothing in the cube. If we denote the union of all these cubes 
WQ (for weak cubes) we find that 

/ [V]X^'< Cvh^-'^^ (15) 

JwQ 

for Cv > a constant depending only on V (in particular on the support of V^). 

In each cube where V^in > V^|| W||oo we have < (1 + "\/3)Knin- Each of these cubes we 
fill with the maximal number of disjoint balls of radius 

R = K,hm'-''^'v-:^^'^-\ 

where ^2 > is a constant which we will choose below depending only on e. Note that 

R < K~'-''/'h'/^+'/''{/3hy-''/^VV\\^'+''/' < Vh/2 (16) 

if h is small enough depending only on V and e (recall that is bounded from above in 
terms of V). In particular, we can then fit at least one ball in the cube. 

In each of these balls we choose a magnetic field according to Proposition 13.31 with 5 
chosen such that 35(1 + 5)'^ = e. Let B : ^ M? he the sum of all these disjointly 
supported magnetic fields and let A be a corresponding vector potential. Let i^'^^^and i?j for 
i = 1,2,... denote the balls (there are only finitely many) and their radii. We can according 
to Proposition 13.31 for each i find a normalized ipi supported in B^^^ such that 

T ■ (-zW + A)^Pi\^ - [ IV^.I V + P [ B^< Csh^l3'~'R;'-' - K.in,i 

J JB(i) 

— ^nin,i/2 
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if K2 is large enough depending on 6 , i.e., on e. Here l^in,i is the minimum of V in the cube 
containing the ball B'^'^\ 

Since all the constructed balls are disjoint we have that P = Yli is an orthogonal 

projection and hence 

E{(3, h, V) < Tr [{{a ■ {-ihV + A)f - V) P] + (3 J 

i i 

< -2Ceh-\l3h)-''+^' [ V^-', 



for some constant > depending only on e. We have here used that in each of the cubes 
the balls take up a certain fraction bounded below of the volume and that V < {1 + \/3)Knin- 
Let us emphasize that K2 was chosen depending only on e. This allows us to ensure that f|T6|) 
is satisfied if h is small enough depending on V and e. Finally, from f[T5|) we can choose h so 
small depending on V and e that 

[v]t^ < \ I [vt-\ 

This proves the claim f[M|) . □ 



4 Semiclassics for weak fields: Proof of Theorem 2.2 



Proof of Theorem We discuss only the Pauli case, the Schrodinger case is similar but 
much easier and is left to the reader. In the remaining part of the proof we will omit the 
superscript P. 

For the upper bound on h^E{(3, h, V) we choose A = in the definition ([1]) and then we 
have E{I3, h, V) < E{oo, h, V) and the second equality in ([8]) is just the usual non-magnetic 
semiclassical asymptotics. 

For the lower bound, we first remark V > can be assumed and that it is sufficient to 
prove the result for V G C^(M^) by a standard approximation argument. The error between 
V E L^/^ nL^ and its C^- approximation V can be made arbitrarily small in the || ■ II5/2 + || ■ lU 
norm. Thus the replacement of V with V can be controlled by using the magnetic Lieb- 
Thirring inequality ffTTj) and borrowing a small part of the kinetic energy and the magnetic 
energy. For more details, see Section 5.4 of |ES3] (with the only modification that instead of 
(5.58) of [ESS] use (^). 

Secondly, as discussed in Appendix R] we may in ([T]) replace J^^ |V x by f^^ |V ® Ap, 
where the last integrand contains all derivatives. 
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Thirdly, we may replace E{I3, h, V) with a localized version of the total energy. Let < 
< 1 be a smooth function with supp0* C B{1) and 0* = 1 on 5(1/2), where B{r) 
denotes the ball of radius r centered at the origin. Denote by (prix) = (j)*{x/r). Using a 
partition of unity, 0^ + 77^ = 1, and the IMS localization, we have 

Tr [n{A) - V]^ > Tr [MTh{A) -V - /i^/,)^,]^ + Tr [r/.(T,,(yl) -V- h^Ir)Vr]_, (17) 

where := (V^,.)^ + (V?7,.)^ is supported in the shell {r/2 < \x\ < r}. The second term is 
bounded by the magnetic Lieb-Thirring inequality similarly to f[T^ . More precisely, for any 
e > there is a sufficiently large r = > 1 such that 

Tr [vr{Th{A) -V- h'lr)vr]_ > Tr (T,(A) - V). 

>-Ch'^ [ V""'^ -Ch-\ph)-''( f V^)-^f\VxA\ 

J\x\>r ^ J\x\>r ' 2 y 



I / IV X A| 



2 



holds ii h < h^, where V{x) := V{x) + /i^/^. Here we used the integrability conditions on V 
and that Ph is bounded from below. 

It is therefore sufficient to give a lower bound on the first term in ( 1T71) . more precisely, we 
have 

E{f3, h, V) > -eh~^ + inf [Tr [(l){Th{A) - W)(f\ _ + ^ /" | V ® Af\ , (19) 

A L 2 J^i J 

where we set W := V + h?Ir and = 0^ for brevity. 

To estimate the right hand side of ( IT9|) . we will follow the argument of Section 5 of |ES3] . 
We choose a length L with h < L < h^^^. Let il^ := B(r + L) be the L- neighborhood of 
Q := B{r). Let Qk := {?/ G : \\y - k\\oo < L/2} with k G (LZ)^ fl Ql denote a non- 
overlapping covering of B{r) with boxes of size L. In this section the index k will always 
run over the set (LZ)'^ fl Ql. Let be a partition of unity, "^^^l = 1, subordinated to the 
collection of boxes Qk, such that 

supp C (2Q)fc, \V^k\<CL-\ 

where {2Q)k denotes the cube of side-length 2L with center k. Let ^k be a cutoff function such 
that = 1 on {2Q)k (i.e. on the support of ^k), supp^^ C Qk '■= {3Q)k and |V^fc| < CL~^. 

Let (y4)fc = IQfcl"^ Jq^ a, Ak := (A— (A)fc)^fc and Bk := Vx74fc, then by Poincare inequality 
we have 

' Bl< [ |V®A,f <c/ |V®Ap + CL-2/ 

>^Qfc jQk JQk 

<C [ \V®A\\ (20) 



13 



From the IMS localization with a phase function satisfying hVipk = (^fc) we have 
Tr [(j){n{A) - W)<j)] _ + ^ / 52 = inf Tr (i<P[mA) - W]<p) + /3 / | V ® 



>inf ^^(^) 



(21) 



fcG(LZ)3nt7i 



with 



+ C0/3 / |V®Ap 



with some universal constant cq and after reallocating the localization error. In fl2T|) the 
infimum is taken over all density matrices < 7 < 1. We also reallocated the second integral 
to account for the finite overlap of the cubes Q^. We introduce the notation 

^fc := Co/3 / |V ® Af. 
JQk 

Let [H\q denote the operator H with Dirichlet boundary conditions on the box Q. For 
each fixed box Qk we apply the magnetic Lieb-Thirring inequality [LLS] together with fl20l) 
to obtain that for any density matrix 7 



£fc(7)>Tr [Th{Ak)-W -Ch'L 



2 r-2i 



>-Ch-' [ [W + Ch^L-'f'-c( f [W + Ch^L-^y^Uh-^ [ Bi\"\:F, 
jQk ^JQk ^ ^ JQk ^ 



IV® + J- 



> -Ch~'L 



3 r3 



using h < L and I3h — t- 00. The constants C depend on 

Let S C [LTjY n VLl denote the set of those k indices such that 

J-fc < Cih-^L^. (22) 

holds with some large constant Ci. In particular, by choosing Ci sufficiently large, we have 

£k{,l) > 0, for all /c ^ 5 and for any 7. (23) 

We use the Schwarz inequality in the form 

Th{A - {A)u) > -(1 - ek)h'A - Ce-\A - (A),.)^ 
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with some < e^,. < |, then 



+ Tr 



l.[-e,h'A-Ce-,\A-{A),f]l. 



(24) 



We will show at the end of the section that 



Tr 



;k[-{^-2ek)h'A-W -Ch'L-% 



> Tr 



Ch-\eu + h'L-')\Q 



(25) 



Using ([23D and dlSD, 

inf^£fc(7)>inf5^^fc(7) 

^ k kes 

> 5^ Tr [<j)^k[~h^A - W]^k<l>] ^ + 

k k&S 

>^infTr h,^,^,(f)[-h'A-W]<p\+Y,1^k 

k 

>Tr [0(-/i2A_iy)0]_ + J^p, 



(26) 



kes 



kes 



with 



V, := Tr 



[skh'A - C6^\A - {A)kyU - Ch-'lQklisk + h^L-^) + Tk 



(27) 



In the last step in 0261) we used that for any collection of density matrices 7^, the density 
matrix J2k ^klk^k is admissible in the variational principle 



Tr [(/.(-/i2A-1^)0]_ = inf |Tr7[-/i2A-M/] : < 7 < 1, |. 
We estimate for G 5" as follows 

Vk>-Ce^'h-^ [ (A - {A)kf - Ch-'\Qk\iek + h'^L-'') + 
JQk 

> - Cel%-^(5~'>l'L^I^:Fl" - Ch~^\Qk\{ek + /^'^-'). 



(2^ 



(29) 



In the first step we used Lieb-Thirring inequality, in the second step Holder and Sobolev 
inequalities in the form 



/ {A-{A)kf<CL''^( [ |V®A|2) 



5/2 
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We choose 

and using the a priori bound (122|) . we see that 

Ek < Ch^^L{l3hy^'^. 

Thus, choosing 

L = h{l3hY'^\ (30) 
we get Ek < C{(3h)~'^/^ < 1/3 as (3h ^ oo. With this choice of Ek and L we have from fl29l) 

^^fe > J^fe - C{hl3)-^/^J^l'^ - Ch-''L\(3h)-^"' > -Ch-^L\(3h)-^'\ (31) 
Summing up flMl) for all /c and using that 

we obtain from fl2B]) and flM]) 

inf y £fc(7) > Tr ^(-/i^ A - W)(f)] - Ch^^/Shy^/^r^ 
2 



> - 



157r2 



(l + o(l)) W'/^ -Ch-''{f3h)-^/V 



(32) 



using the standard semiclassical asymptotics for Tr [0(— /i^A — Vr)0]_ > Tr (— /i^A — iy)_. 
Together with (12T]) this proves the required lower bound for the second term in ( !T9|) . The 
difference between J VT^/^ and J V^^'^ is negligible as /i — > 0. Letting first /i — > together 
with — 7- oo and then letting 5 — ?■ we obtain the lower bound in (jS]). 

Finally, we prove f l25|) . Let 7 be a trial density matrix for the left hand side of fl25|) . We 
can assume that 



> Tr 



l<P^k[-{l - 2EkWA -W- Ch^L-^]ik 



Recalling that Ek have 



>Tr 
>Tr 



+ Tr 



^cPik[-\h''A-W-Ch^L-^ 



^<P^k[--h^A + l\ik<p\ -Ch-^ i [W + l + Ch-'L-^f^ 



(33) 



where we used Lieb-Thirring inequality. Thus, using h < L, we have 

1 



Tr 



70a[-g^'A + l]a0j <Ch-'\Qk\ 
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with a constant depending on W. Therefore 



Tr 



70a[-(l - 2ek)h^A -W- Ch^L^^j^k 



> Tr 



Now (123]) follows by the variational principle. □ 



(34) 



A Equivalent forms of the energy 

We will consider the equivalence of the total energy where we have different restrictions on 
the vector potentials. We allow the energy to possibly have an extra localization. So we end 
up considering 



where ip & C 



E^{A)=TT[^{Th{A)-V)ij]. + ^ f \VxA\ 
satisfies < ip < 1. Similarly, we define 

E^{A) = Tr [i^iThiA) - V)^;]. + /3 / | V ® A| 



where |V ® = Yl^j=i l^i^iP- Some natural domains of definition are given below: 



AeL^ 



: J\WxA\^< oo}, 



2 = C^iW',W'), 



V4 = {AeVi : V ■ A = 0}. 



In the case of E^, the expression J |V x in "Di should be replaced by J |V ® 

We will only 
be interpreted as 



(35) 



(36) 



We will only assume that V G Lj'Q^(M^). The trace in the above expressions should then 



N 



(37) 



where {^jj^^ runs over all orthonormal subsets of C^(M^). If this infimum is different from 
— oo, it implies in particular that the quadratic form of ip{Th{A) — V)ip defined on is 
semibounded from below. In that case (1571) will be equal to the trace of the negative part of 
the Friedrichs extension of this quadratic form, thereby justifying the notation. 
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Proposition A.l. We have for alii, j G {1,2,3,4}, 



inf E^{A) = inf E^{A). (38) 

Notice though that we do not prove that one can both impose compact support and zero 
divergence at the same time. 

Proof. Consider first the Ey^. Clearly, V2 C V^, C Vi (using the Sobolev inequality to get 
the last inclusion) which implies corresponding inequalities for the energies. We will now 
prove that iniAi^Vi Ex{A) > inf^gp^ -^x (^)- But for any A E Vi and any finite collection 
{4>j} C C^(]R^) we can get arbitrarily close to 

N 



5^(</.,#(T,(A) - V)^^J<PJ) + P f\VxA\' 
.7=1 



by simultaneously approximating A in L^-norm and V x A in L^, by a vector field. 
Therefore 

inf Ey(A) = inf E^iA) = inf Ey,(A). 

AgVi A&V2 A£Vs 

Clearly voIa^v^ (A) < inf AgD4 E^ (A). We will prove that M^eVi Ex [A) > ^AeVi E^ (A), 
thereby establishing the equality for all four energies E^. Let A G C^(]R^,M'^) and {4)j}jLi C 
C^(M3). Then B = V x A e L'^{R^) and therefore there exists A' G with V x A' = B 
(see [FLL] ). It follows that there exists 77 with A — A' = Vt] and therefore (since Arj = V ■ A), 
r] G C°°(M3). But then 

N N 

i=i i=i 

which establishes the desired inequality. Since J|V®Ap = / |Vx/lp + J|V-Ap, the same 
arguments give the identities for the E,^ versions of the energies. 

Finally we prove that iniA^Vi Ex{A) = inf^gx'4 E^{A). But this is obvious since the field 
energies are identical when V ■ A = 0. □ 



B Self-generated magnetic fields lower the energy 

In this appendix we will show that self-generated magnetic fields may indeed decrease the 
energy, i.e., inequality ([3]). 

Proof of For the Pauli operator we already remarked this fact as a consequence 
of Theorem 12.11 Alternatively, it also follows from the instability in f HTj) in Appendix O 
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below, since the non-magnetic Hydrogen atom is stable. For the Schrodinger operator this 
statement was essentially proved in |ELV] (see also |FLWj ) by considering the perturbative 
regime as a small magnetic field is turned on. A simple first order perturbation argument shows 
that the lowest eigenvalue increases quadratically in B. In a spherical geometry the higher 
non- magnetic eigenvalues are degenerate and some of them carry non-trivial current. These 
eigenvalues will split linearly when a small magnetic field is turned on. To see this explicitly 
we can consider a spherically symmetric harmonic oscillator in a constant magnetic field, i.e., 
V{x) = |xp and A{x,y,z) = {By/2, —Bx/2,0) with B > constant. The eigenvalues of the 
operator (— iV — A)^ + |xp are (see |Fo] ) 

e(ni, ns, ng) = (ni + ria + l)Vl + B^ + (ng + 1/2) + (rii - n2)B 

with ni,n2,n3 G N. Thus as an explicit example 

Tr ((-iV - + |xp - 5/2)_ = 3Vl + B^ - 4 - 5 

which of course explicitly decreases as a small B is increased from zero. It is now clear that 
we can find A G C^(M^; M^) which approximates A such that 

Tr ((-zV - Af + |a;p - 5/2)_ < Tr (-A + |xp - 5/2)_ = -1 

and hence for /3 > sufficiently small 

Tr ((-iV - Ay + - 5/2)_ + (3 J |V x Ij^ < Tr [-A + - 5/2]_. 

□ 

C Stability conditions 

Using the argument in jFLL] it is easy to show the following stability result on the one-electron 
energy. 

Proposition C.l. Let V G L^M^) with V+ G L^R^) n L^/'^{R^). Then for all ip G H^{R?) 
with ||^||l2 = 1 we have 

{^lj,{Tl{A)-V)^l^)+P Jb'>0 

if (]~^h~'^\\V+\\s and h~'^\\V+\\3/2 are (universally) small enough. 
Proof Let C5 > be the Sobolev constant, i.e., / \V'ip\'^ > Cs\\i>\\l. Since 

T^(A) = {-ihV + Af + h(T-B 
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we estimate for all < £ < 1 

>cshh\mi-m~'e'h' j iv^r 

We will also use that for p > 3/2 we have the Holder inequality 

/ v\i^\'<\\v4\mT\mV''- 

We consider two cases. 

Case 1: < 2CsP- We set e = 1 and p = 3/2 above and find since 



2 



{i;,{T^{A)-V)^) + {3 1 B'>\C,K^ 



^-Il^+ll3/2||^||^ 



from which it follows that the energy is non- negative if ||V+||3/2 < j'^- 
Case 2: ||^||6 > 2Cs/3. Let e = 2Cs^\\^\\q^ < 1 and = 3. Then 

Hence the energy is non-negative if ||V^+||3 < C'^h'^p. □ 

It follows immediately from this proposition that the one-electron energy E^{l3,h,V) is 
finite if V+ G L^(M^). In fact, all we have to argue is that [V — e] + ||3 and /;,~^|| [V^ — e] + ||3/2 

can be made small enough by choosing e > large enough. In this way — e can be made a 
lower bound on Eq. Since Vj^ G we can of course make || — e] + II3 arbitrarily small. Using 
that [V - efl^ < (2e-i)3/2[y - e/2]^ we can do the same with \\[V — e] + ||3/2. 

This stability criterion is essentially sharp. In fact, applying the method of proof as in the 
proposition above and the construction of zero-modes in |LY] it was proved in |FLL] that for 
the Coulomb potential V{x) = c|x|~^ there is a critical value 7cr such that the one-electron 
energy satisfies 

E^{(3, h, V) > -00 if 7cr/3/i^ > c (39) 

and 

E^{(3, h, V) = -00 if 7cr/3/i' < c. (40) 

Since E^ < Eq it is clear that (l40l) implies that even for the cutoff Coulomb potential V = 
[c\x\~^ — 1]+ we have 

E^{(3, h, V) = -00 if 7cr/3/i^ < c. (41) 
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However there is also a value 7^^, > such that 

E^{/3, h, V) > -00 if y^^f3h^ > c (42) 

This stability statement follows, e.g., by localizing in an appropriate ball and then follow 
the proof of Lemma 2.1 jES3j for the inner regime (with the choice of Z = h""^, 5 = Z^^^, 
D = RZ^/^ and levra^ = (5-^). In the outer regime the operator has a compactly supported 
bounded potential (that includes the localization error) so its energy is controlled by the 
magnetic Lieb-Thirring inequality as in f[T2|) . 
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